Neutron EDM from Electric and Chromoelectric Dipole Moments of Quarks 
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Using QCD sum rules, we calculate the electric dipole moment of the neutron d n induced by 
all CP violating operators up to dimension five. We find that the chromoelectric dipole moments 
of quarks di, including that of the strange quark, provide significant contributions comparable in 
magnitude to those induced by the quark electric dipole moments di. When the theta term is 
removed via the Peccei-Quinn symmetry, the strange quark contribution is also suppressed and 

d n = (1 ±0.5) |l.le(d d + 0.5d u ) + 1.4(d d - 0.25d u )}. 



Experimental discovery of CP violation in flavor- 
conserving channels would provide a clear indication of 
new physics around the electroweak scale. This stim- 
ulates continuing experimental efforts to detect electric 
dipole moments in elementary particles and heavy atoms 
[1]. The extraordinary precision, d n < 6 x 10 -26 ecm 
[2], obtained in measurements of the electric dipole mo- 
ment (EDM) of the neutron, allows us to probe energy 
scales inaccessible in direct collider experiments. The in- 
formation on fundamental CP violating parameters, such 
as CP-odd phases and masses of superpartners in super- 
symmetric (SUSY) models, is encoded in the coefficients 
of CP-odd effective operators 0W, evolved down to 1 
GeV. However, our ability to extract this information 
depends critically on the quality of calculations of the 
neutron EDM induced by these operators. 

In this Letter, we report the results of the first sys- 
tematic analysis of d n , induced by all CP violating oper- 
ators with dimension 4 and 5, within the QCD sum rule 
framework [3]. The complete set {O^'} of operators at 
the 1 GeV scale includes the theta term and the electric 
and chromoelectric dipole moments of quarks (EDMs and 
CEDMs): 

SC = - V q(m q + ie gl5 )q + 6 G ^GG 
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where F nv and G M „ are the electromagnetic and gluonic 
field strength tensors. With these sources, the neutron 
EDM can be written as a linear combination of the coef- 
ficients in Eq. (1), 

d n = d n (9) + df M (d u ,d d ,d s ) + d^EDM fa J d ^y (2) 

The first term d n {6) arises from the dimension four 9— 
term in (1), where 9 = 9 q + 9q is the physical combi- 
nation. In the absence of e.g. the Peccei-Quinn (PQ) 
compensation mechanism [4], this is the most impor- 
tant source term, and the calculation of d n (9) within 



QCD sum rules was performed in [5, 6]. Here we 
shall present the calculation of the dimension five con- 
tributions, d^ DM (ci 9 ) and d^ EDM (J g ), which are required 
for phenomenological analyses. In fact, the results of 
d n ]UM (d q ) can be extracted from calculations of the ten- 
sor charges of light quarks in the proton [7] . At the same 
time, earlier calculations of d n rEBM (d q ) in QCD sum rules 
[8] predict an unexpected suppression of CEDM contri- 
butions, contrasting with chiral loop estimates which in- 
dicate sizable effects [9]. Therefore, the main focus of 
this work is the crucial task of calculating the contribu- 
tions of CEDMs, which also depend on the presence or 
absence of the PQ mechanism [10]. 

The starting point for the calculation is the correlator 
of currents r) n (%) with quantum numbers of the neutron 
in a background with nonzero CP-odd sources and an 
electromagnetic field F uu , 



n(Q 2 ) = 



= i 



d 4 



xe^ x {0\T{ Vn (x)r] n (0)}\0)cp.F, (3) 



where Q 2 — —p 2 , with p the current momentum. 

In the presence of CP violating sources it is necessary 
to take into account mixing between the neutron current 
and its CP conjugates. Thus we parametrize the inter- 
polating current r\ n in the form, 



Vn = (Ji ± Ph) ± ieqp(»i + Ph) 



(4) 



where the two conventional neutron interpolators, j± = 
2£abc(daC~/ 5 u b )d c and j 2 = 2e ahc (d^Cu b )^ b d c , are com- 
bined with their CP conjugates, i\ = 2e a b c (d a r Cub)d c and 
i2 = 2e a bc(d a r j5Cub)j5d c . The parameters e and j3 in (4) 
play rather different roles. Firstly, e reflects mixing in- 
duced by the CP violating source and will be calculated 
below. Thus to linear order in the sources, 

(VnVn) = 0'J)f,cp + ieqp(3 i + i 3)F +0(e c p) 2 , (5) 
where, 



2 0'iji - Mi) 



(6) 
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The second parameter (3 reflects the existence of the 
two interpolators ji and j'2, and would disappear from 
the result of an exact calculation. It can therefore be 
used to optimize the convergence of the operator product 
expansion (OPE) of (3). Within the sum rules formalism, 
one has the imperative of suppressing the contribution of 
excited states and higher dimensional operators in the 
OPE, and thus its convenient to choose (5 to this end. 
We shall therefore keep (3 arbitrary, and optimize once 
we have knowledge of the structure of the sum rule. 

In analyzing the correlator (3) an additional consider- 
ation is that when CP-symmetry is broken by a generic 
quark-gluon CP-violating source, the coupling between 
the physical neutron state, described by a spinor v, and 
the current -q n acquires an additional phase factor, 



\'h 



\N) = Ae M75/2 w, 



(7) 



where the coupling naturally decomposes as A = Ai +/3A 2 
which in part motivates our parametrization above, with 
the same parameter /3 for the CP conjugate currents i\ 
and i 2 . With regard to electromagnetic form factors, 
the unphysical phase a in (7) can mix electric (d) and 
magnetic (/i) dipole moment structures and complicate 
the extraction of d from the sum rule. As discussed in 
[5] , there is a unique tensor structure contributing to the 
neutron double pole term which is not contaminated by 
this phase, namely {Faj5,]/>}. Therefore, it is this struc- 
ture that we shall use to construct the EDM sum rule 
from (3). An additional advantage of using this struc- 
ture comes from the absence of unknown vacuum tensor 
polarizabilities [7] in the calculation of d^ r>M (d q ). 

We now proceed to study the OPE associated with (3). 
To next-to-leading order in the OPE the relevant classes 
of diagrams we need to consider are shown in Fig. 1 ((a), 
(b) and (c)). Diagrams of the form (d), although suffering 
no loop factor suppression, are nonetheless suppressed 
due to combinatorial factors and the small numerical size 
of ({qq)) 2 relative to the natural scales in the problem. 
(This partially explains why the results of previous anal- 
yses [8] indicated a suppression of d^ EI}M (d q )). 

The vacuum structure is conveniently encoded in a 
generalized propagator expanded in the background field 
and the associated condensates. The condensates are 
then parametrized in terms of various susceptibilities de- 
fined as [11]: 

(q<r^q) F = XqF^uiqq); 9(q(G a ^t a )q) F = KgF^fiq) 
g(qGaq) = -m 2 ; 2g(q l5 (G%t a )q) F = i^F^qq), (8) 

and henceforth we follow [11] and assume that \q = X e q 
etc., with flavor independent susceptibilities x, K >£,- The 
dependence on CP violating parameters is either explicit, 
e.g. 9 q and d q , or implicit in certain vacuum condensates 
in the case of d q and will be discussed shortly. 

The relevant contributions to the OPE are exhibited 
in Fig. l(a-c), but the explicit expressions are quite un- 
wieldy, and we shall defer full details [12], and simply 




(a) 



(b) 





(d) 



FIG. 1: Various classes of diagrams contributing to the CP-odd 
structure {Fa-ys,^)}. 



present the resulting OPE expression. In momentum 
space we find 



n(0 2 ) = -%£-te){F<n*,t} 



64tt 2 



+ 



16tt 2 p 2 



In 



■^)(qq){F<n5,t}[*l + *i], (9) 
V / 



corresponding to the first two nontrivial orders in the 
OPE. The relevant contributions from 9 and the CEDMs 
are contained in, 

ttW = 4(1 + Pfxam d P d - (1 + f3) 2 Xu m u P u 

+2(1 -f3 2 )m*( X u + Xd)(Pu-Pd), (10) 



and 



\d d [a+(3 + 2p + 3p 2 )-a+(l-P 2 ) 

-a-(l - pf] + l -d u [2a+(l - P 2 ) - a+(l + P?] , (11) 

where we have defined P q = 9 q — i ((775(7) qp / (qq) , to* = 
to u to < j/(to u + TO<j), and a ± — 2k ± £ for concision. In 
obtaining this result it is crucial to take into account the 
mixing with the CP conjugate currents in (5). The pa- 
rameters controlling the mixing are the unphysical phases 
9 a — 9 q and the quark mass difference m d — m u . Here, as 
expected, the addition of the mixing terms ensures the 
absence of the unphysical phases 9q — 9 q . 
The quark EDM contributions are 



n q = d d [10 + 6/3 2 ] - d u [3 + 2/3 - P 2 



(12) 



At subleading order, terms with a logarithmic momen- 
tum dependence arise and involve an infrared cutoff [Iir. 
The relevant contributions are, 



(1 - (3) 2 m d e d P d 
+ {l-(3 2 )m«{e u + e d ){P d -P u ) (13) 

1 



12 



(1 - l3) 2 ml(d d - e d d d ) 



(14) 
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The next problem to address is the calculation of 
the vacuum matrix elements P u and Pd in Eq. (9). 
These terms require the evaluation of correlators of the 
form j d 4: y(q"f^q(x),iSC(y)), where 5£, given in Eq. (1), 
involves in particular the 0-term and the color EDM 
sources which may be extracted from the vacuum at lead- 
ing order in the background electromagnetic field. The 
case of 9 was discussed at length in [5, 13, 14]. Here we 
shall concentrate on the CEDMs and evaluate these cor- 
relators in chiral perturbation theory, saturating them by 
7To and rjs exchange, assuming the decoupling of a heavy 
singlet state, and incorporating the effect of the chiral 
anomaly [15]. We then obtain, 



m u (d)Pu(d) 



m*ml I d u (d) - d d ( u ) 



m d ( u ) 




where we have neglected terms of 0(d( u ^/m s ), assum- 
ing an approximate proportionality of the CEDMs to the 
quark masses, i.e. dd/d s ~ md/m s <C 1. 

With these results in hand, we now turn to optimiza- 
tion of the OPE expression (9). Recall that there are 
generally two motivated approaches for fixing the mixing 
parameter (3: (1) at a local extremum; or (2) to minimize 
the effects of the continuum and higher dimensional op- 
erators. Since we are restricted here to only the first two 
orders in the OPE, we make use of the second option, as 
in [5], and set j3 = 1 to cancel the subleading infrared 
logarithmic terms, which are ambiguous due to the cut- 
off. This procedure mimics the original motivation for 
(3 = — 1 in the CP-even case [16]. 

On the phenomenological side of the sum rule we have 



n phon = -{^ 7 5,. 



X 2 d n m n 
(P 2 - ml ) 2 



A 



mt 



.(16) 



We retain here the double and single pole contribu- 
tions, the latter corresponding to transitions between the 
neutron and excited states, but the exponentially sup- 
pressed continuum contribution will be ignored. In (16) 
A = Ai +/3A 2 and A is an effective constant parametrizing 
the single pole contributions. A more detailed analysis 
will be presented elsewhere [12], and we have verified that 
the inclusion of a continuum does not affect the analysis 
below. 

After a Borel transform of (9) and (16), and using (3 = 
1 as discussed earlier, we obtain the sum rule, 



\ 2 m n d n + AM 2 



M 4 
'32^ 



X/M' 



Ax) 



'(3=1+^(3=1 



+ 7T 



/3=1 



(qq) 

+ 0{M 2 ) 1 (17) 



where the contributions now take the elegant form, 



Ax) 
T /3=i 



4 [4x d m d Pd - Xum u P u ] 



(£,«) _ 



0=1 



'/3=1 



4d d a^ - d u a+ 



4 [4d d - d u ] . 



(18) 

(19) 
(20) 
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FIG. 2: The neutron EDM function u(M 2 (GeV 2 )) is plotted ac- 
cording to the sum rules (a) and (b). 



It is remarkable that the contribution of the u quark is 
— 1/4 that of the d quark, which is precisely the combi- 
nation suggested by the SU(6) quark model! 

In order to analyze the sum rule (17), we first deter- 
mine the coupling A using the sum rules for the tensor 
structures 1 and j> in the CP even sector (see e.g. [17] 
for a recent review). Following [5], we construct two sum 
rules. Firstly, (a): we extract a numerical value for A 
via a direct analysis of the CP even sum rules. This 
analysis has been discussed before and will not be repro- 
duced here (see e.g. [17]). One uses (3 = —1, and obtains 
(27r) 4 A ~ 1.05 ±0.1. As an alternative, (b): we extract A 
explicitly as a function of (3 from the CP-even sum rule 
for and substitute the result into (17) choosing (3=1. 

The conventional approach which we shall adopt here 
is to assume that A is independent of M, and thus the 
left hand side of (17) is linear in M 2 provided that A 
is constant in the appropriate region for the Borel pa- 
rameter. For case (b), the latter point may be verified 
explicitly. The function v(M 2 ), given by 



v{M 2 ) = 



d n + 



AM' 



A 2 



ni. 



(21) 



is then determined by the right hand side of (17). 

The two sum rules described above for v a and Vb are 
plotted in Fig. 2. v(M 2 ) is to be interpreted as a tan- 
gent to the curves in Fig. 2. For numerical calcula- 
tion we make use of the following parameter values: For 
the quark condensate, we take (qq) = —(0.225 GeV) 3 , 
while for the condensate susceptibilities, we have the val- 
ues ml = 0.8 GeV 2 [18], X = -5.7 ± 0.6 GeV" 2 [18], 
£ = -0.74 ±0.2 [8], and k = -0.34 ±0.1 [18]. 

One observes that both sum rules have consistent ex- 
trema suggesting that our procedure for fixing the param- 
eter (3 is appropriate. Furthermore, the differing behavior 
away from the extrema implies that for consistency we 
must assume A to be small. One then finds d n as given 
by v{M 2 ~ 0.55GeV 2 ). This low scale is characteristic of 
CP violating effects, and convergence of the OPE is ap- 
parently not in danger due to the combinatoric suppres- 
sion factors which arise in extracting the CP violating 



4 



sources from higher dimension operators. 

Extracting a numerical estimate for v{M 2 ~ 
0.55GeV 2 ) ~ 0.043 from Fig. 2, and determining an ap- 
proximate error due to higher order corrections and the 
sum rule analysis, we find the result [20]: 



d n = (0.4 ±0.2) 



X(4e d - e„) ( m*6> - \mld s — 
2 m s 



, 1 2/7 7 Ae d m d 
+ 2* m o(«<i - d u ) 



e,.m, 



m u + m d 



+ o( 4 ^ a d - duQ-t,) + ( 4d d - d u ) 



(22) 



for the neutron EDM. The CEDM contributions are sig- 
nificant and comparable in magnitude in fact to the ef- 
fects induced by quark EDMs. 

For most phenomcnological applications, it is nec- 
essary to invoke a PQ symmetry, which removes the 
dominant contribution d n {6). However, while setting 
d n (0) = 0, this symmetry induces additional CP vi- 
olating terms through linear contributions to the ax- 
ion potential [10]. In particular, the axion poten- 
tial has the form V ~ -6 2 K - 26K', where K = 
i(a s /8Tr) 2 (GG,GG) is the topological susceptibility, and 
K' = i(a s /8ir)(GG, <5£ CEDM ) are correlators arising from 
the CEDM sources. This linear shift in the axion poten- 
tial then leads to an "induced" 9-term with coefficient 
(sec e.g. [15]) 



K 1 



_ _-n_ _ m \ - 
ind ~ K ~ 2 ^ 



2 



d a 



q—u,d,s 



m a 



(23) 



which importantly is independent of any of the specific 
details of the axion mechanism. 

Including these contributions, we observe the complete 
cancellation of the term proportional to the strange quark 
CEDM, and (22) takes a remarkably simple form 



d p n Q = (0.4 ±0.2) 



Ad d - d u + \p(m\{^e d d d - e u d u ) 



+ -{M d aj - d u a+) 



(24) 



Substituting numerical values for the condensates, we ob- 
tain 



= (1±0.5) 



1(95)1 



(225McV) 3 
l.le(d d + O.bdu) + lA(d d - 0.25d u ) 



(25) 



The result for d EUM (d q ) agrees well with lattice calcula- 
tions of quark tensor charges of the proton [19], and with 
the naive quark model estimate. Note that an overall 
factor of (qq) combines with the light quark masses from 
short-distance expressions for di and di to give a result 



~ f%m 2 (1 + 0(mr Utd )) thus reducing the uncertainty due 
to poor knowledge of the quark masses and condensates. 
Further progress in calculating d n would require more 
elaborate analysis of the sum rules in order to reduce the 
error in the overall coefficient, while the relative coeffi- 
cients of different terms in the square brackets of (25) are 
likely to remain the same. 

In conclusion, we have presented the first systematic 
study of the neutron EDM induced by all CP violat- 
ing sources of dimension four and five within QCD sum 
rules. We observe that the contributions from the CEDM 
sources, including that of the strange quark, are large - 
in fact as large as those associated with the quark EDM 
sources. This result has significant implications for the 
analysis of the EDM constraints imposed on the super- 
symmetric and other models of CP violation. 
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